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I.  INTRODUCTION 


Chemical  lasers  offer  a unique  means  of  obtaining  high 
intensity,  coherent  radiation.  Unlike  other  types  of  lasers 
which  require  large  amounts  of  external  power,  the  chemical 
laser  obtains  lasing  action  through  the  energy  released  in  a 
chemical  reaction.  Conventional  chemical  lasers  employ  either 
H2  and  P2  or  D2  and  F2  as  the  principal  reactants  in  the 
system.  Normally,  the  lasing  action  occurs  in  the  mixing 
region  between  a stream  of  gas  containing  primarily  H2  or  D2 
and  another  stream  containing  F2  or  F and  a diluent,  e.g.. 

He.  The  performance  of  the  laser  depends  on  the  ability  of 
the  configuration  to  mix  and  react  the  gases  over  a distance 
sufficiently  large  that  optical  energy  may  be  efficiently 
extracted.  Thus,  the  proper  mixing  and  subsequent  chemical 
reactions  producing  excited  chemical  species  is  one  of  the 
most  important  problems  in  the  design  of  a good  chemical  laser. 

Current  chemical  lasers  operate  at  low  cavity  pressures, 
5-10  torr,  in  order  to  reduce  the  rate  of  collisional  deacti- 
vation of  the  excited  species.  At  these  low  pressures,  the 
mixing  of  the  reactants  is  predominantly  through  molecular 
diffusion,  a rather  slow  process.  If  the  mixing  could  be 
achieved  at  a more  rapid  rate  then,  conceivably,  the  rate  of 
excitation  would  increase  considerably,  providing  higher  laser 
output  power  and  overall  efficiency.  In  addition,  if  the 
mixing  could  occur  at  a faster  rate,  then  it  may  be  possible 
to  operate  the  laser  cavity  at  higher  pressures  thereby  reducing 
the  problem  of  bringing  the  gas  up  to  atmospheric  pressure 
(supersonic  diffusion). 

A more  rapid  mixing  of  the  reactants  may  be  achieved 
through  turbulence  in  the  flowfield.  The  possible  role  of 
turbulence  in  chemical  lasers  can  best  be  understood  by  exam- 
ining the  structure  of  flames.  In  a laminar  flame  the  reactions 
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between  the  fuel  and  oxidizer  are  confined  to  a very  narrow 
region.  In  this  case,  the  primary  mode  of  mixing  is  through 
molecular  diffusion.  Employing  molecular  diffusion  coefficients, 
the  so-called  flame  velocity  may  be  theoretically  calculated  to 
a reasonable  degree  of  accuracy.  As  the  flow  velocity  is  in- 
creased and  the  flow  becomes  turbulent,  two  observations  can 
be  made.  First,  the  flame  velocity  is  much  higher  and,  secondly, 
the  flame  zone  becomes  very  wide.  The  widening  of  the  flame 
zone  is  an  indication  that  the  dominant  mass  transfer  is  through 
turbulent  mixing.  In  this  situation,  the  width  of  the  flame 
zone  approaches  the  scale  of  the  turbulence.  Unlike  the  laminar 
flame  velocity,  the  prediction  of  the  turbulent  flame  velocity 
is  very  difficult  because  of  the  lack  of  a sufficient  under- 
standing of  turbulent  reacting  flows. 


The  chemical  laser  presents  a somewhat  similar  situation. 
Laminar  mixing  produces  a very  thin  zone  of  lasing  action  be- 
tween the  streams  of  reactants.  As  the  flowfield  becomes  tur- 
bulent, the  lasing  zone  becomes  wider.  Conventional  chemical 
lasers  operate  in  the  transition  regime  in  which  the  flow  cannot 
be  characterized  as  fully  laminar  or  fully  turbulent.  The 
prediction  of  turbulence  in  the  transition  regime  requires 
sophisticated  techniques.  Turbulent  mixing  models  which  do  not 
consider  the  history  of  the  turbulence,  such  as  the  eddy  vis- 
cosity approach  and  the  turbulent  kinetic  energy  approach,  may 
lead  to  considerable  errors  in  predicting  laser  performance. 


This  situation  is  particularly  important  in  flow  situations  in 
which  turbulence  is  introduced  through  various  boundary  layer 
trip  devices.  In  these  situations  It  is  important  that  a 
theoretical  approach  be  capable  of  treating  variable  turbulent 
scales,  intensities,  and  distributions  of  energy  among  the  u , 
v , and  w components  of  velocity. 


Although  turbulence  may  enhance  the  mixing  of  the  reactants, 
it  may  have  adverse  effects  on  the  rates  of  the  individual 


chemical  reactions.  In  the  turbulent  flame  example  cited 
above,  the  use  of  a laminar  chemistry  model  leads  to  consider- 
able error  in  the  prediction  of  the  flame  width.  The  laminar 
chemistry  model  means  that  the  rate  of  a particular  reaction 
is  determined  only  by  the  average  concentrations  of  the  re- 
acting species  and  the  local  average  temperature.  Although 
the  local  average  species  concentrations  have  been  determined 
through  turbulent  mixing,  the  reactants  may,  in  fact,  not  be 
mixed  on  a molecular  level  which  is  required  for  the  species 
to  react.  This  may  be  illustrated  by  considering  the  initial 
mixing  of  black  and  white  paint.  After  a few  swirls  of  the 
mixing  paddle,  we  can  define  an  average  concentration  of  the 
two  paints  over  a small  region.  But,  clearly  these  two  are 
not  mixed  on  the  molecular  level.  It  is  this  unmixedness 
which  is  of  primary  importance  in  treating  the  effects  of 
turbulence  on  chemical  reactions. 

Fully  realizing  the  problems  associated  with  the  theoret- 
ical prediction  of  turbulent  reacting  flowfields,  A.R.A.P. 
began  developing  a turbulence  code  which  could  provide  the 
types  of  information  which  are  required  to  determine  the  effects 
of  turbulence  on  the  rates  of  individual  chemical  reactions. 

In  addition,  it  was  important  that  the  turbulence  model  be 
invariant  since  it  may  be  employed  in  flow  situations  for  which 
no  turbulence  measurements  are  available  to  calibrate  the  model. 
A. R. A. P.'s  second-order  closure  technique  is  the  only  turbulence 
model  which  meets  these  necessary  criteria.  This  report  will 
briefly  describe  the  A.R.A.P.  Reacting  Shear  Layer  (RSL)  code. 

In  addition,  we  will  present  numerous  computations  showing  the 
effects  of  initial  turbulence  scale,  intensity,  and  distribution 
of  turbulent  energy  on  the  generation  of  excited  chemical 
species  in  a chemical  laser.  Furthermore,  we  will  indicate  the 
error  which  can  be  introduced  when  the  adverse  effects  of  the 
turbulence,  the  unmixedness,  is  ignored. 
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II.  ANALYSIS  OP  TURBULENT  REACTING  FLOWS  - 
DESCRIPTION  OF  THE  RSL  CODE 

A complete  second-order  closure  program  for  the  analysis 
of  chemically  reacting  turbulent  flowfields  is  being  developed 
at  A.R.A.P.  The  philosophy  of  our  second-order  closure 
approach  and  details  of  the  technique  have  been  previously 
described  in  a number  of  reports  and  publications  (Refs.  1,2). 
For  the  sake  of  completeness  a detailed  description  of  the 
models  and  the  computer  program  have  been  included  in  this 
report . 

The  theoretical  analysis  of  turbulent  flows  is  made 
difficult  by  the  closure  problem;  namely,  that  the  number  of 
unknowns  in  an  exact  time-averaged  formulation  always  exceeds 
the  number  of  equations.  The  problem  becomes  more  complex 
in  the  case  of  turbulent  reacting  flows  because  of  the  larger 
number  of  fluctuating  quantities  that  are  involved  and  the 
appearance  of  certain  new  and  difficult  scalar  correlations. 
One  commonly  used  approach  for  the  solution  of  the  closure 
problem  requires  the  development  of  models  for  the  extra 
unknown  correlations  to  obtain  a closed  set  of  equations  which 
can  then  be  solved  by  a variety  of  methods.  One  simple 
modeling  approach  that  has  long  been  used  in  the  study  of 
turbulent  flows  is  an  eddy  viscosity  or  first-order  closure 
approach.  Models  are  developed  for  second-order  correlations , 
e.g.,  uiuj  > u|T’  , etc.,  in  terms  of  the  mean  quantities 
of  the  flow.  This  approach  completely  ignores  the  dynamics  of 
the  turbulence  and  considers  it  to  be  a "local"  quantity;  that 
is,  the  turbulence  at  a point  in  the  flow  is  completely  deter- 
mined by  the  local  values  of  the  mean  flow  variables.  The 
procedure  has  been  quite  successful  in  a number  of  engineering 
problems,  but  it  suffers  from  the  need  to  determine  the 
empirical  constants  for  each  class  of  flows  from  direct  exper- 
imental measurements.  It  appears  unlikely  that  an  "invariant" 
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or  "universal"  eddy  transport  theory  can  be  developed. 

Further,  there  are  a number  of  turbulent  flows  of  interest 
for  which  the  dynamics  (or  history)  of  the  turbulence  is 
important  in  determining  the  nature  of  the  turbulence  at  a 
point  in  the  flow.  In  this  respect,  the  first-order  closure 
or  eddy  transport  concept  is  not  a sufficiently  powerful 
tool  with  which  to  study  the  motion.  Some  examples  are  flow 
problems  involving  large-scale  atmospheric  motions,  flow 
involving  chemical  reactions,  and  the  developing  regions  of 
classically  self-similar  flows.  In  the  past,  most  engineering 
calculations  of  reacting  flows  have  been  made  by  ignoring  the 
interaction  between  the  turbulence  and  chemistry.  However, 
it  is  generally  accepted  that  the  interaction  between  turbu- 
lence and  chemical  reactions  is  of  considerable  importance 
in  many  practical  flow  problems,  and  first-order  closure 
methods  are  not  capable  of  handling  such  nonequilibrium  flow 
regions.  Second-order  closure  techniques  are  required  to 
analyze  these  flowfields. 

Considerable  advances  have  been  made  in  recent  years  in 
the  development  of  second-order  closure  models  for  nonreacting 
turbulent  flows  which  have  demonstrated  improved  predictive 
capability  compared  to  the  eddy  viscosity  or  first-order 
closure  models.  A.R.A.P.  has  been  one  of  the  leading  groups 
in  the  development  of  multi-equation,  complete  second-order 
closure  modeling  of  turbulent  flows.  The  application  of  these 
higher-order  closure  procedures  to  turbulent  reacting  flows 
is  a natural  extension  of  current  technology.  In  a second- 
order  closure  procedure  one  develops  transport  equations  for 
second-order  correlations  and  solves  them  along  with  the 
equations  for  the  mean  fluid  dynamic  variables  subject  to  the 
appropriate  initial  and  boundary  conditions.  These  equations 
contain  third-order  and  higher-order  correlations  which  are 
modeled  in  terms  of  the  lower-order  correlations.  The  research 
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effort  of  a number  of  groups  is  directed  toward  establishing 
suitable  models. 


Considerable  success  has  been  achieved  in  the  prediction 
of  incompressible  and  compressible  nonreacting  flows  using 
second-order  closure.  The  formulation  of  the  problem  differs 
from  group  to  group.  Launder  and  Spalding  (Ref.  3),  Bradshaw 
(Ref.  4),  and  Saffman  (Ref.  5)  consider  only  equations  for 
the  turbulent  kinetic  energy  and  turbulence  scale  in  addition 
to  the  mean  flow  equations.  Donaldson,  Varma,  and  others  at 
A.R.A.P.  have  developed  a multi-equation  closure  scheme  (Refs.  6, 
7,  8)  which  retains  individual  equations  for  all  the  important 
second-order  correlations  in  the  problem  and  provides  much  more 
detailed  information  about  the  turbulence  field  than  the  simpler 
schemes.  A. R. A. P.'s  approach  has  been  to  select  the  models  in 
the  simplest  possible  way  consistent  with  tensor  invariance  and 
available  experimental  data.  The  model  must  remain  invariant 
with  respect  to  changes  in  the  flow  geometry.  The  invariant 
modeling  is  an  important  feature  in  that  the  objective  is  to 
obtain  a model  applicable  to  a wide  class  of  fluid  flows.  Data 
from  simple,  extensively  studied  flows  can  then  be  used  to 
obtain  the  general  model  parameters  for  use  in  more  complex 
flowfields  where  experimental  measurements  are  difficult. 

This  model,  with  the  use  of  the  same  set  of  model  parameters, 
has  demonstrated  excellent  agreement  with  experimental  data 
for  a wide  class  of  flowfields,  such  as  free  shear  layers 
(Ref.  6),  flat  plate  boundary  layers  (Ref.  6),  axisymmetric 
jets  and  wakes,  the  planetary  boundary  layer  (Ref.  7),  and 
the  decay  of  turbulent  vortices  (Ref.  8).  Finson  (Ref.  9) 
has  also  used  multi-equation  models  with  good  results. 

The  use  of  second-order  closure  modeling  for  reacting 
flows  is  still  in  its  infancy.  In  addition  to  the  normal 
difficulties  associated  with  the  second-order  closure  analyses 
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of  turbulent  flows,  the  presence  of  finite  rate  chemical 
reactions  introduces  the  problem  of  the  proper  treatment  of 
the  combined  effects  of  concentration  and  temperature  fluc- 
tuations. At  A.R.A.P.,  Donaldson  and  Hilst  (Ref.  10)  studied 
low  heat  release  chemical  reactions  and  demonstrated  the 
importance  of  keeping  track  of  the  species  correlations. 

Order  of  magnitude  errors  can  be  made  in  the  prediction  of  the 
reaction  rate  in  practical  problems  if  classical  chemical 
source  expressions  (based. on  mean  quantities  alone)  are  used 
and  the  species  and  temperature  correlations  ignored.  Spalding 
(Ref.  11)  formulated  an  "eddy  break  up"  model  valid  for  large 
turbulence  Reynolds  numbers  and  for  very  fast  chemical  reactions 
to  model  the  effect  of  turbulence  on  chemical  reactions.  Chem- 
ical kinetic  effects  are  ignored  and  the  rate  of  the  reaction 
is  taken  to  be  proportional  to  dissipation  rate  of  the  macro- 
scale eddies.  Recently,  Spalding  (Ref.  12)  has  proposed  a 
LaGrangian  approach  to  the  modeling  of  turbulent  combustion 
which  attempts  t p edict  the  development  of  the  probability 
density  functions  in  the  flow. 

The  major  new  problem  in  the  analysis  of  compressible 
reacting  flows  is  the  modeling  of  the  large  number  of  higher- 
order  correlations  of  scalar  fluctuations  arising  in  the  chem- 
ical source  terms.  The  general  approach  being  used  by  many  of 
the  research  groups  active  in  the  field  is  to  model  the  proba- 
bility density  function  (pdf)  for  the  scalars.  These  include 
Rhodes  et  al.  (Ref.  13),  Bray  and  Moss  (Ref.  1*0,  Libby  (Ref  15), 
and  Borghi  (Ref.  16). 

Most  of  these  approaches  still  make  many  simplifying 
assumptions  to  make  the  problem  tractable.  Bray  and  Libby  use 
the  fast  chemistry  assumption.  The  approach  of  Borghi  is  very 
similar  to  ours,  but  he  replaces  the  joint  pdf  of  species  and 
temperature  by  a one-dimensional  pdf  which  is  a very  restrictive 
assumption.  The  A.R.A.P.  approach  models  the  joint  pdf  of  all 
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the  scalars  In  turbulent  reacting  flows.  A procedure  has  been 
developed  for  constructing  the  pdf  using  the  available  Infor- 
mation In  a second-order  closure  analysis.  The  model  is  called 
the  "typical  eddy"  model  and  involves  representing  the  pdf  by 
a set  of  delta  functions  of  variable  strengths  and  positions  in 
the  scalar  phase  space.  The  strengths  and  positions  of  these 
delta  functions  at  every  point  in  the  flow  are  determined  from 
the  predicted  values  of  the  means  and  second-order  correlations 
at  that  point.  Kewley  (Ref.  17)  has  recently  used  our  concept 
of  the  "typical  eddy"  model  in  a study  of  turbulence  effects  in 
a chemical  laser.  The  details  of  the  models  and  the  RSL  program 
code  for  the  analysis  of  general  multi-step  exothermic  chemical 
reactions  in  turbulent  shear  flows  are  discussed  below. 

Basic  Equations 

The  general  conservation  equations  for  a turbulent,  com- 
pressible, multi-component  reacting  flow  system  are  presented 
in  Appendix  A.  The  following,  commonly  used  assumptions  have 
been  adopted  for  the  derivation  of  these  equations: 

1.  The  gas  is  a continuum. 

2.  No  body  forces  are  present. 

3.  Radiant  heat  transfer  is  neglected. 

4.  There  are  no  overall  mass  sources. 

5.  Mass  diffusion  occurs  by  a concentration  gradient 
only.  Fick's  diffusion  law  is  valid.  Thermal  and 
pressure  diffusion  are  neglected. 

6.  Each  component  of  the  gas  mixture  is  thermally  perfect. 

These  equations  contain  a number  of  third-order  and  higher- 
order  correlations  that  have  to  be  modeled  in  terms  of  the 
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means  and  second-order  correlations  to  obtain  a closed  set  of 
equations . 


Second-Order  Closure  Models 


Modeling  of  Pressure  Fluctuation  Correlations 


Tendency-toward-isotropy  terms 
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In  these  models,  A is  the  turbulent  scale  length  to  be 
determined.  The  major  change  compared  to  the  models  used  in 
References  6 and  19  is  the  addition  of  terms  involving  the 
mean  strain  9um/9xn  . Reference  19  suggests  P2  = 0.1  . 
Constants  ?1  and  P3  are  still  under  investigation. 


Models  for  Other  Correlat ions 


Velocity  diffusion  terms 

Following  previous  successful  A.R.A.P.  studies  (Refs.  6, 
19),  we  adopt  a gradient  diffusion  model s as  follows: 
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In  Eqs.  (6)  and  (7),  <t>  ’ represents  any  scalar  variable 


( T1 


etc.  ).  Vf,  was  established  equal 


> > (x  5 H »c 

to  0.1  by  comparing  model  results  with  experimental  data  in 

free  jets  and  shear  layers.  Due  to  the  lack  of  detailed 
experimental  data  in  compressible  and  reacting  flows,  the  same 
diffusion  constant  is  used  as  an  initial  estimate  for  velocity 
diffusion  of  all  correlations. 


Dissipation  terms 


Lewellen  et  al.  (Ref.  20)  suggest  that  isotropic  dissipa- 
tion is  more  suitable  than  the  anisotropic  dissipation  model 
used  by  Donaldson  (Ref.  19). 
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X is  the  turbulent  microscale  length 
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Detailed  studies  of  the  flat  plate  incompressible  boundary 
layer  (Ref.  21)  have  led  to  a slight  change  in  a : 


a = 3.25 
b = 0.125 
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Models  for  Scalar  Correlations 


A.R.A.P.'s  "Typical  Eddy”  Model 

The  major  new  difficulty  in  second-order  closure  computa- 
tions of  turbulent  reacting  flows  is  the  modeling  of  scalar 
correlations,  such  as  kfa’  , k ' a * 6 ' , a ' a * g ' , etc.,  that 
appear  in  the  equations.  Recently,  we  have  developed  a "typical 
eddy"  model  that  is  very  promising  for  the  modeling  of  all  such 
scalar  correlations.  Details  of  the  development  of  the  model 
are  given  in  References  22  and  23.  A brief  description  of  the 
model  is  given  below. 

The  construction  of  the  "typical  eddy"  model  is  an  attempt 
to  define  the  joint  probability  density  function  (pdf)  for  all 
the  scalars,  using  the  available  information  from  a second-order 
closure  analysis.  The  pdf  is  represented  by  a set  of  delta 
functions  of  variable  strengths  and  positions  in  the  scalar 
phase  space.  The  location  of  some  of  the  delta  functions  is 
decided  on  physical  grounds.  The  positions  of  the  other  delta 
functions  and  the  strengths  of  all  the  delta  functions  at  every 
point  in  the  flow  are  determined  from  the  values  of  the  means 
and  second-order  correlations  obtained  at  that  point  from  the 
solution  of  their  transport  equations.  The  degrees  of  freedom 
in  the  model  are  limited  by  the  number  of  independent  parameters 
available  in  the  second-order  closure  analysis.  Our  complete 
second-order  closure  program  for  a three  species  flow  system 
provides  information  for  13  independent  first  and  second-order 
correlations  — 5 , 6 , a ' 6 ' , a ' y ' , B ' y ' , a ’ p ’ , B ' P ' , 
p » p ' , h , h * h * , a * h * , g * h*  , and  p ’ h ' — so  the  complete 
model  can  have  no  more  than  13  parameters.  Figure  1 shows  the 
model  that  we  have  selected  after  extensive  studies.  The  model 
has  exactly  13  parameters  to  be  determined.  The  typical  eddy  is 
considered  to  have  the  following  structure.  For  a fraction  e-j.  » 
the  eddy  contains  only  the  species  a . The  cell  has  a sensible 
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enthalpy  ha  and,  since  for  a single  species  flow  from  a con- 
stant enthalpy  source  the  model  must  allow  enthalpy  fluctuations, 
h ' h ' , due  to  boundary  conditions  and/or  viscous  heating,  we 
add  a Ah  for  half  of  the  cell.  For  a single  species  flow, 
the  enthalpy  structure  has  just  two  parameters  ( hQ  and  Ah  ) 
to  be  calculated  from  the  values  of  h and  h * h ' obtained 
from  the  transport  equations.  In  the  model  currently  being  used 
in  our  reacting  shear  layer  programs,  the  same  Ah  is  used  for 
all  the  cells.  For  a fraction  £2  > the  eddy  contains  only 
species  6 with  an  enthalpy  hg  , and  for  a fraction  e^  » 
the  eddy  has  only  the  pure  species  y and  has  the  corresponding 
enthalpy  h^  . For  fractions  e^  , , and  eg  , the  eddy 

is  assumed  to  contain  a and  6 , a and  y , and  g and  y , 
respectively,  in  a state  of  molecular  mixedness.  There  appear 
to  be  two  equally  reasonable  choices  for  the  proportions  of  the 
species  in  these  cells.  One  can  consider  the  species  to  be 
present  in  either  equal  amounts  or  in  amounts  proportional  to 
their  average  values  at  that  point  in  the  flow.  The  current 
program  retains  both  of  these  options  at  the  present  time,  and 
further  test  runs  will  indicate  the  more  desirable  model.  Cor- 
responding average  values  are  used  for  the  enthalpies  of  the 
cells.  Finally,  an  e 7 fraction  of  the  eddy  contains  all 
three  species  a , 8 , and  y in  a molecular  mixture  of 

proportions  <1  » <2  > and  1 - - <2  • The  enthalpy  of 

cell  7 is  h 7 . The  model  has  13  unknowns  ( e-^  through  eg  , 

<1  , <2  > ha  , hg  , hY  , hy  , and  Ah)  (note  X c±  = l) 

to  be  determined  from  the  13  available  independent  correlations 
by  matching  the  moments  of  the  model  to  the  values  predicted 
by  the  equations.  Thus,  the  model  can  be  constructed  at  any 
desired  time  from  the  second-order  closure  turbulence  calcu- 
lation in  progress.  Once  the  species  and  enthalpy  distribution 
functions  have  been  established,  the  corresponding  distributions 
for  the  temperature,  density,  and  the  reaction  rate  can  be  con- 
structed. The  complete  nonlinear  Arrhenius  rate  expression  can 
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be  used  for  k ; there  is  no  need  to  expand  the  exponential 
term.  All  scalar  correlations  required  to  close  the  system 
of  second-order  correlation  equations  can  now  be  calculated 
from  this  joint  probability  density  function. 

The  construction  of  the  complete  "typical  eddy"  model 
described  above  requires  the  solution  of  a set  of  non linear 
algebraic  equations.  A simplified  "typical  eddy"  model  was 
also  proposed  (Refs.  22,  23).  This  model  does  not  consider 
the  density  fluctuation  correlations  in  the  construction  of 
the  model  and,  therefore,  can  only  have  9 independent  parameters. 
The  simplified  model  results  in  a linear  set  of  equations  for 
the  determination  of  the  cell  sizes  in  the  species  probability 
density  function,  and  analytical  closed  form  solutions  have 
been  obtained.  In  contrast,  one  has  to  resort  to  numerical 
solutions  in  the  case  of  the  complete  model.  The  simplified 
"typical  eddy"  model  was  easier  to  implement  for  the  initial 
testing  of  the  concept  of  the  model. 

The  simplification  involves  the  neglect  of  the  density 
correlations  in  the  setting  up  of  the  model  and,  therefore, 
should  be  valid  for  flows  where  density  changes  are  small.  The 
simplified  "typical  eddy"  model  has  been  shown  (Ref.  1)  to  be 
satisfactory  for  low  heat  release  reacting  flows.  However, 
just  how  far  one  could  use  this  simplified  scheme  could  only 
be  tested  by  actual  computations  and  comparison  to  experimental 
measurements.  Our  studies  on  chemical  laser  flows  and  hydrogen- 
air  flames  where  significant  density  changes  are  involved 
suggest  that  the  more  complete  "typical  eddy"  model  has  to  be 
used,  as  was  expected,  for  these  flows. 

An  even  simpler  model  has  been  used  In  the  calculations 
for  the  DF  chemical  laser  reported  here.  In  this  model  all 
third-order  and  higher-order  scalar  correlations  are  set  to 
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zero.  This  is  designated  as  the  "secondary"  model,  as  distin- 
guished from  the  primary  model  which  is  the  "typical  eddy" 
model.  Second-order  correlations  such  as  T' a ? , k ’ a * , and 
others  for  which  we  do  not  carry  transport  equations,  are  now 
obtained  up  to  second-order  accuracy  by  expansion  of  explicit 
thermodynamic  relations  between  the  instantaneous  variables. 

When  the  "typical  eddy"  model  is  used,  these  conversion  rela- 
tions are  not  necessary  as  the  pdf  for  p,T  and  k can  be 
constructed  directly  from  the  pdf  of  oij_  and  h using  the 
thermodynamic  relationships.  The  "secondary"  model  was  initially 
used  during  the  model  development  stage  but  some  recent  calcu- 
lations (Ref.  1)  have  indicated  that,  at  least  for  some  simple 
reacting  flows,  the  two  models  may  predict  quite  similar  results. 
It  can  be  shown  quite  easily  that  the  "secondary"  model  is  not 
correct  in  a number  of  limiting  cases.  For  example,  in  the 
reaction  end  limit  for  a system  of  three  species  undergoing  the 
one  step  reaction  a + 8 -*•  y , 

aB  = 0 
a8  + a ' 8 ’ = 0 

oB*  + a'8  + a ' 8 ' - a ' 8 ' = 0 
a a ' 8 ' + a'a'  6 + a ' a ' 8 ' = 0 

a ' a ' 8 ' = -a  a'8'  -8  a’a’ 

■ • 5!(>  - W) 

t 0 in  general 

Thus,  setting  the  third-order  scalar  correlations  to  zero 
is  not  a proper  model  and  the  agreement  between  the  results  from 
the  "secondary"  model  and  the  "typical  eddy"  model  must  be  simply 
coincidental.  The  computations  for  the  DF  chemical  laser  will 
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be  repeated  later  with  the  "typical  eddy"  model  and  the  com- 
parison of  the  two  sets  of  results  should  be  quite  interesting. 
It  is  an  interesting  question  why  the  results  for  some  flows 
that  have  been  studied  are  insensitive  to  the  choice  of  the 
model  for  the  higher-order  scalar  correlations. 

Multi-Step  Chemistry  Procedure 

The  reacting  shear  layer  program  and  the  "typical  eddy" 
model  were  originally  developed  for  use  with  three  species 
( a , 0 , and  y ) flowfields,  and  only  permitted  the  single- 

step  reaction  a + 0 -*•  y . The  concept  of  such  an  overall 
reaction  may  be  entirely  sufficient  in  some  combustion  problems, 
but  for  most  problems  it  is  necessary  to  use  a set  of  elementary 
reaction  steps  involving  many  species  to  correctly  describe  the 
system  and  to  obtain  results  in  better  agreement  with  experi- 
ments. This  conclusion  is  particularly  valid  with  respect  to 
trace  species.  Further,  the  elementary  reactions  have  been 
studied  extensively  and  the  reaction  rate  parameters  are  better 
known.  Ideally,  one  would  like  to  keep  track  of  each  species 
and  the  correlation  between  each  species  and  the  other  species 
as  well  as  the  correlation  with  certain  flow  properties.  Unfor- 
tunately, this  desire  introduces  an  enormous  number  of  equations 
which  must  be  solved  in  an  implicit  fashion.  The  resulting 
computer  time  would  be  prohibitive  to  the  point  that  the  code 
would  simply  not  be  used  for  engineering  calculations.  There- 
fore, it  was  decided  to  construct  a system  for  handling  a 
finite  number  of  species  and  reactions  within  the  present  frame- 
work. 


To  understand  the  basic  philosophy  of  our  reaction  scheme, 
it  is  useful  to  consider  the  expression  for  the  rate  of  change 
of  a species,  a , due  to  chemistry  in  a situation  in  which 
turbulence  may  be  present.  It  can  be  shown  that  the  expression 
reduces  to: 
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= -kcig  -[ka ' g ' + akUS7  + 3k_roTr  + k'a'  8'  ] 

Laminar  Turbulent  Contribution 

In  the  absence  of  turbulence,  the  rate  of  change  of  the  species 
is  simply  given  by  the  laminar  expression.  When  turbulence  is 
present,  all  of  the  other  terms  must  be  considered.  The  com- 
putation of  these  terms,  and  the  other  necessary  correlations, 
is  the  reason  the  computer  time  becomes  excessive.  In  the 
laminar  case,  such  correlations  are  not  required.  Thus,  the 
extension  of  laminar  flow  situations  to  hundreds  of  chemical 
species  and  reactions  is  straightforward  since  it  only  requires 
the  addition  of  another  equation  for  the  conservation  of  mass 
of  each  added  species.  In  the  turbulent  case,  the  addition  of 
a new  species  requires  the  addition  of  at  least  three  equations 
for  the  correlations  with  the  primary  flow  variables  plus  an 
equation  for  the  correlation  with  each  of  the  other  species  in 
the  flow.  From  a practical  viewpoint,  and  the  desire  to  main- 
tain reasonable  computation  times,  it  was  decided  to  "track" 
only  three  species:  a , 6 , and  y . 

In  any  diffusive  reacting  flow,  one  usually  has  two  in- 
itially separate  streams  of  fuel  (a)  and  oxidizer  (g).  These 
streams  are  now  considered  to  be  known  mixtures  of  a number  of 
chemical  species.  The  streams  mix  and  react  following  a large 
number  of  elementary  reactions  and  form  a mixture  of  product 
species  designated  as  y . We  assume  that  (1)  the  composition 
of  mixtures  a and  g is  fixed  for  the  entire  flowfield. 

There  are  no  internal  reactions  among  the  species  in  these  two 
reactant  mixtures.  (2)  All  the  product  species  that  compose 
y are  molecularly  mixed.  There  are  internal  reactions  within 
the  y species  mixture.  The  composition  of  the  mixture  y 
varies  at  different  points  across  the  flowfield. 
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Consider  the  following  two-reaction  system  as  an  illus- 
tration of  the  procedure.  The  complete  set  of  reactions  that 
have  been  used  in  our  studies  of  the  DF  chemical  laser  are 
presented  in  Table  1. 

^1 

D2  + F — i-  DF  + D 
^2 

F2  + D 4-  DF  + F 
a D2  known 

B He  + F + F2  known  proportions 
Y He  + D2  + F2  + D + F + DF 

The  interaction  between  turbulence  and  chemistry  only  has 
to  be  taken  into  account  for  reactions  between  a and  6 , 
a and  y , and  6 and  y . Species  within  y are  assumed 
to  be  molecularly  mixed  and,  therefore,  species  correlations 
do  not  have  to  be  considered  for  the  internal  reactions.  For 
convenience,  we  rewrite  the  reaction  system  as  given  below, 
tagging  each  chemical  species  with  the  mixture  that  it  is  a 
part  of: 

D2(a)  + F(B)  - DF(y ) + D(y) 

D2  (a)  + F(y)  -*•  DF(y ) + D(y) 

D2(y)  + F(8)  - DF(y)  + D(y) 

D2(y)  + F(y)  -*■  DF(y)  + D(y)  internal  reaction 

F2(B)  + D(y)  ■*>  DF(y ) + F(y) 

F2(y)  + D(y)  DF(y)  + F(y)  internal  reaction 
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The  chemical  source  terms  for  the  component  species  can 
now  be  written.  For  example. 


dD2(a) 


dt 


= - k„ 


D2(a)F(B)  + D£(a)F'(B) 


- k. 


D2(a)F(y)  + D2 ( a ) F ’ ( Y ) 


dD2(y) 


dt 


= - k2  D2(y)F(B)  + ( y )F ' ( B ) 


- kx  D2(y)F(y) 


, etc  . 


Transport  equations  are  solved  for  the  second-order 

fTT7 


correlations  a * 0 ' 


a*Y' 


and 


The  additional 


important  assumption  is  made  that  quantities  like  D£(a)F'(0) 


D£(a)F'  (y) 


etc.  can  be  simply  calculated  from  these  corre- 


lations as  being  proportional  to  the  local  mean  composition 
of  the  mixture.  The  validity  of  this  assumption  has  to  be 
examined  further,  but  it  should  not  lead  to  large  errors.  The 
rate  at  which  a or  0 react  to  form  the  product  species  y 
is  governed  by  the  rate  of  the  fastest  reaction.  In  the  above 
example,  0 is  being  lost  due  to  the  consumption  of  both 
F ( 0 ) and  F2(0)  • The  reaction  that  causes  the  greater  loss 
of  6 is  the  controlling  step.  If  at  some  point  in  the  flow, 
the  reactions  of  F ( 0 ) lead  to  the  larger  loss  of  0 , then 
to  maintain  the  constant  composition  of  the  0 species  mixture, 
corresponding  amounts  of  F2 ( 0 ) and  He (6)  now  become  part 
of  the  y species.  The  composiiton  of  y has  to  be  calculated 
at  each  point  across  the  flowfield.  The  procedure  operates 
within  the  framework  of  three  overall  species  mixtures,  and 
all  the  models  and  equations  developed  earlier  can  be  used.  The 
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only  addition  to  the  program  is  the  solution  of  the  mean 
species  conservation  equations  for  all  the  elemental  species. 
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Numerical  Solution  of  Equations 

The  use  of  the  models  for  various  third-  and  higher-order 
correlations,  as  described  earlier,  in  the  set  of  equations  for 
the  means  and  second-order  correlations  results  in  a closed  set 
of  equations.  Equations  for  the  mean  variables  p , u^  , h , 
a , and  6 and  the  second-order  correlations  u|u j , h* h'  , 

u|h'  , a ' 8 ' , , a'Y*  , B ' y ' , up*7  , u |8 ' , u^p ' , h'aT  , 

and  h * B ' are  solved  in  the  present  program.  The  use  of  the 
basic  shear  layer  assumptions  leads  to  a set  of  23  independent 
coupled  parabolic  partial  differential  equations.  The  computer 
program  actually  solves  a total  of  30  equations;  the  redundant 
equations  provide  a check  on  the  mass  conservation  in  the  program 
and  on  the  accuracy  of  the  numerical  scheme.  The  numerical 
integration  of  the  equations  is  performed  by  a forward-time- 
centered-space , quasi-implicit , upwind,  finite  differencing 
scheme.  The  nonlinear  terms  are  handled  by  quasi-linearization, 
that  is,  by  evaluating  a portion  of  these  terms  at  the  known 

position,  leaving  only  a linear  term  containing  one  of  the 

unknowns.  The  linearized  finite-difference  equations  are  solved 
by  the  general  tridiagonal  algorithm;  but,  Instead  of  solving 
a single  large  matrix  (which  would  be  very  time  consuming),  the 
equations  are  grouped  into  smaller  matrices  and  the  system  is 
solved  in  10  separate  passes.  The  program  is  running  satis- 
factorily on  our  Digital  Scientific  Corporation  META-4  computer. 

The  program  has  the  capability  of  handling  both  fixed  and 
free  shear  layer  flows  for  planar  and  axisymmetric  geometries. 

The  initial  profiles  can  be  provided  using  a card  or  tape  input 
or  the  program  can  generate  appropriate  smooth  initial  profiles 
given  the  properties  of  the  two  streams.  If  the  turbulence 
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properties  at  the  initial  station  are  not  known , we  typically 
input  a "spot"  of  turbulence  at  the  initial  station  to  start 
the  turbulent  calculations.  The  program  uses  a fourth-order 
polynomial  expression  to  calculate  the  temperature  dependence 
of  the  specific  heat  of  various  species.  The  molecular  tranS' 
port  properties  y , y*  , k , and  D are  calculated  using 
a power  law  expression. 


III.  MODEL  AND  PROGRAM  VERIFICATION 


The  principal  objective  of  developing  a higher-order 
closure  turbulence  model  than  the  first-order  closure  eddy- 
viscosity  models  is  to  obtain  a more  "universal"  model;  a 
model  which  can  be  used  for  analysis  of  data  for  a wide  class 
of  flowfields  with  the  use  of  the  same  invariant  set  of  model 
constants.  Such  a model  can  then  be  used  with  a greater  de- 
gree of  confidence  for  predictive  calculations  of  flow 
problems  for  which  experimental  data  are  not  available  or  are 
difficult  to  obtain.  Second-order  closure  models  have  shown 
significant  promise  of  being  such  a predictive  tool. 

The  fluid  mechanical  turbulence  models  being  used  in  the 
A.R.A.P.  second-order  closure  program  have  already  been  ex- 
tensively tested  by  comparison  of  program  predictions  to  experi- 
mental data  in  a variety  of  basic  flow  geometries.  Previous 
publications  have  discussed  studies  on  flat  plate  boundary 
layers  (Ref.  6), the  planetary  boundary  layer  (Ref.  7),  and 
axisymmetric  jets  and  wakes  (Ref.  24).  Model  and  program  veri- 
fication studies  that  are  of  interest  in  connection  with  the 
analysis  of  the  planar  mixing  flow  in  DF  chemical  lasers  are 
our  analyses  of  a two-dimensional  wake  and  a heated  planar  jet 
(Ref.  1)  and  the  mixing  of  two  different  species  in  a shear 
layer. 

The  experimental  measurements  in  the  wake  of  a thin  plate 
were  made  by  Chevray  and  Kovasznay  (Ref.  25).  Detailed  measure- 
ments of  the  mean  axial  velocity  profile  and  the  Reynolds  stress 
were  made  using  hot  wire  probes.  The  experimental  measurements 
are  designated  as  Test  Case  14  in  data  collected  at  the  1972 
NASA  Free  Turbulent  Shear  Flows  Conference  (Ref.  26).  Figure  2 
compares  the  theoretical  RSL  program  predictions  for  the  mean 
axial  velocity  profiles  at  three  axial  stations  with  the 
measurements.  The  computations  are  generally  in  good  agreement 
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with  the  data  over  the  entire  profile.  There  is  some  dis- 
crepancy between  the  predictions  and  the  experiments  as  to 
the  centerline  velocity  at  x = 20  cm  , the  reason  for  which 
is  not  known  at  this  time.  Figure  3 shows  the  predicted  pro- 
file for  the  Reynolds  stress  correlation  u ' v'  at  x = 50  cm  . 

At  this  axial  station,  the  computations  are  in  good  agreement 
with  the  data.  Further  downstream,  at  x = 240  cm  , the  peak 
value  of  the  computed  u'v'  is  about  25%  lower  than  the  experi- 
mental data,  but  the  profile  shape  is  in  good  agreement. 

Davies,  et  al.  (Ref.  27)  and  Uberoi  and  Singh  (Ref.  28) 
have  carried  out  extensive  measurements  of  mean  and  turbulence 
quantities  in  heated  plane  turbulent  jets.  Figure  4 presents 
our  program  predictions  for  the  decay  of  the  mean  velocity  and 
temperature  on  the  jet  centerline  and  compares  the  predictions 
to  the  data  of  Ref.  27.  The  decay  of  the  mean  velocity  is  in 
excellent  agreement  with  the  data  and  the  linear  graph  shows 
the  self-preserving  nature  of  the  turbulent  jet.  At  larger 
axial  distances,  the  turbulence  falls  off  and  the  decay  rate 
deviates  from  the  straight  line.  The  decay  of  the  centerline 
temperature  is  well-known  to  be  faster  than  the  velocity 
as  shown  here.  The  slope  of  the  calculated  temperature 
(6  = T - Tamb^ent)  decay  rate  is  slightly  smaller  than  Davies’ 
experimental  results.  However,  the  overall  agreement  of  the  pro- 
gram predictions  with  experimental  measurements  is  quite  satis- 
factory. In  the  self-preserving  region,  the  calculated  values  of 
u7uTl/'2/umax  are  0.16  at  the  centerline  and  a peak  value  of 
0.20  , in  agreement  with  Davies'  results  of  0.18  and  0.21  , 

respectively,  as  shown  in  Figure  5.  Figure  6 shows  the  results 

1 /?  _ 

for  the  profile  of  0'0'  /0max  at  x/b  = 20  . Our  program 
predicts  0.22  at  the  center  and  a peak  value  of  0.31  at 
y /£g  = 0.8  . Davies’  measured  values  are  significantly  lower 
and  the  profile  is  virtually  flat.  Uberoi  and  Singh  have  also 
made  temperature  correlation  measurements  in  a plane  heated  jet 
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and  these  are  in  much  better  agreement  with  our  predictions, 
and  are  also  more  reasonable  in  shape  considering  the  profile 
of  the  velocity  fluctuations.  The  discrepancy  between  the 
measurements  of  Davies  and  Uberoi  has  to  be  investigated 
further.  In  both  Figures  5 and  6,  the  predicted  profile  shape 
is  in  good  agreement  with  the  measurements. 

Recent  measurements  in  a two-species  mixing  layer  have 
been  reported  by  Konrad  (Ref.  29).  These  measurements  are 
part  of  a very  elaborate  and  detailed  study  of  this  basic 
flowfield  by  the  group  under  the  direction  of  Prof.  Roshko  at 
Caltech.  Konrad  has  made  measurements  of  the  mixedness  corre- 
lation a ' 8 ' /a8  or  aTar/oi(  1-a)  in  a shear  layer  of  velocity 
ratio  = 0.38  and  consisting  of  two  streams  of  different  species 
but  the  same  density.  Figure  7 shows  a comparison  of  the  RSL 
program  predictions  for  the  flowfield  with  the  experimental 
measurements.  There  is  good  agreement  in  the  central  turbulent 
region  of  the  flowfield.  The  theoretical  predictions  for  the 
edge  regions  of  the  flow  are  not  symmetrical,  which  is  some- 
what surprising  and  is  believed  to  be  due  to  the  difficulty  of 
accurately  calculating  the  term  a ' 8 ' /a 6 when  a or  8 is 
very  small.  This  requires  some  further  study.  The  disagree- 
ment between  the  theoretical  predictions  and  the  experiments 
may  also  be  due  to  the  intermittent  nature  of  the  flow  in  this 
region  and  the  observed  large  structures  in  the  turbulent  flow. 
The  theoretical  models  do  not  include  intermittency  effects  at 
the  present  time. 


All  the  above  computations  and  those  mentioned  earlier  have 
been  carried  out  with  basically  the  same  set  of  model  parameters. 
The  results  demonstrate  that  our  second-order  closure  techniques 
correctly  predict  the  mean  flow  and  details  of  the  turbulence 
correlations  for  a wide  class  of  nonreacting  flows  and  provide 
a level  of  confidence  in  the  models  and  the  programs.  Similar 
tests  have  to  be  carried  out  for  turbulent  reacting  flows  to 
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Figure  7.  Profile  of  the  mixedness  correlation  a ' 6 ' /aB  in  a shear  layer. 


verify  the  models  developed  for  the  scalar  correlations. 
Recently,  computational  results  have  been  presented  for  a 
hydrogen-air  diffusion  flame  (Ref.  30)  that  correctly  predict 
some  of  the  important  features  of  the  flow.  However,  many 
more  such  tests  remain  to  be  made  to  thoroughly  check  the 
models  and  the  programs. 
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IV.  CHEMICAL  LASER  FLOWFIELD  COMPUTATIONS 
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The  A.R.A.P.  turbulent  reacting  shear  layer  (RSL)  program 
has  been  used  to  study  the  mixing  and  chemical  reactions  in 
the  initial  region  of  the  flowfield  of  a supersonic  DF  chemical 
laser.  The  flow  parameters  and  the  initial  geometry  for  the 
studies  are  shown  in  Figure  8.  The  design  data  corresponds  to 
the  MESA  IV  slit  nozzle  built  by  Rocketdyne,  and  was  provided 
by  the  Air  Force  Weapons  Laboratory.  The  flow  parameters  were 
originally  for  a HF  chemical  laser  system.  For  our  studies 
we  simply  replaced  the  H2  by  D2  and  maintained  all  the  other 
flow  parameters  at  the  same  level.  The  nozzle  wall  boundary 
layers  and  the  dead  water  regions  at  the  nozzle  exits  are 
ignored  in  these  studies.  The  analysis  of  the  dead  water  region 
is  particularly  difficult  as  recirculating  flows  are  likely  to 
be  present  and  elliptical  flow  equations  have  to  be  solved  to 
properly  characterize  the  region.  The  effect  of  the  initial 
nozzle  wall  boundary  layers,  however,  can  be  included  in  the 
computations  without  too  much  effort  in  the  future. 

The  Initial  profiles  of  the  mean  variables  are  taken  to  be 
smooth  profiles  of  width  .025  cm  (10%  of  the  total  width)  between 
the  specified  values  on  the  nozzle  centerlines.  The  computa- 
tional region  is  from  the  centerline  of  the  D2  nozzle  to  the 
centerline  of  the  fluorine  + helium  stream.  Suitable  boundary 
conditions  are  used  at  these  symmetry  planes.  The  turbulence 
correlations  u'u'  , v* v»  , w* w*  , and  u* v'  are  specified 
at  the  initial  station  to  begin  the  computations.  Typically, 
a spot  of  turbulence  of  variable  amplitude  and  with  u * u ' * v* v’  = 
w' w'  = 2u ' v*  is  used  at  the  nozzle  exit  in  the  mixing  region 
between  the  two  streams.  The  effects  of  varying  the  initial 
turbulence  amplitude  and  the  distribution  of  the  total  turbu- 
lence kinetic  energy  q2  among  the  three  normal  stresses  has 
been  studied.  The  mixedness  correlation  a ' 6 k / a 5 was  also 
specified  as  nonzero  at  the  initial  station  in  some  of  the  runs. 
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Figure  8.  Flow  configuration 


In  these  runs,  this  correlation  was  set  at  -1.0  across  the 
entire  profile  at  x = 0 ,,  and  then  allowed  to  develop 
according  to  its  transport  equation.  For  initially  unmixed 
streams  of  reactants  a and  B , -1  is  the  proper  value  at 
the  initial  station,  for  a and  6 do  not  coexist  at  any 
point  and  are  thus  perfectly  anti-correlated.  However,  the 
initial  value  of  this  correlation  is  not  of  crucial  importance. 

In  identical  runs  with  a'B'/al  = 0 or  -1  at  the  initial 
station,  the  correlation  reaches  the  same  downstream  value  in 
a very  short  distance,  as  would  be  expected  in  a marching 
solution  procedure.  The  effect  of  the  initial  conditions  is 
wiped  out  in  a distance  of  a few  relaxation  lengths  of  the 
flowfield.  All  the  other  turbulence  correlations  are  set  at 
zero  at  the  initial  station.  The  pressure  is  maintained  con- 
stant at  the  specified  cavity  pressure  throughout  the  flowfield. 

The  primary  objective  of  our  research  program  was  to 
parametrically  study  the  effect  of  various  turbulence  flow 
parameters  — turbulence  kinetic  energy  , turbulence 

distribution  or  isotropy,  turbulence  macro-scale  A , and 
the  turbulence-chemistry  interaction  effects  — on  the  initial 
activation  region  of  the  flowfield  of  a supersonic  DF  chemical 
laser.  Concentrating  on  the  initial  excitation  reactions  for 
the  production  of  excited  DF(v)  species,  a 16-step  reaction 
scheme  was  used  for  the  majority  of  the  test  runs  of  the  program 
reported  in  this  report.  This  reaction  scheme  was  specified 
by  AFWL  and  DARPA  personnel  and  is  listed  in  Table  1.  This  set 
of  reactions  is  only  a small  subset  of  the  complete  reaction 
scheme  as  proposed  by  Cohen  (Ref.  31).  The  excitation  reactions 
up  to  the  fourth  vibrational  level  for  DF  have  been  included. 

The  laser  deactivation  reactions  have  not  been  considered  in 
our  studies.  Toward  the  end  of  our  research  program,  some 
modifications  and  additions  were  made  to  the  reaction  scheme 
listed  in  Table  1 to  evaluate  the  sensitivity  of  the  computations 
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to  the  selected  reaction  scheme.  The  changes  were  to  replace 

the  DFO  in  reactions  5-8  by  2 DF(v)  and  to  add  the 

v=0  -4 

deactivation  reactions  with  helium.  The  modified  reaction  scheme 
used  in  some  of  the  runs  is  listed  in  Table  2.  The  use  of  the 
more  complete  reaction  scheme  led  to  very  minor  changes  in  the 
program  predictions. 

The  results  of  our  studies  are  discussed  in  the  next 
section.  The  computations  using  the  reaction  scheme  of  Table  2 
will  indicate  this  fact.  All  other  runs  have  been  run  using 
the  16  reaction  scheme  detailed  in  Table  1. 
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V.  RESULTS  AND  DISCUSSION 

A majority  of  current  HF  and  DF  chemical  laser  systems 
operate  at  cavity  pressures  of  the  order  of  10  torr.  The  flow 
Reynolds  number  of  these  devices  is  of  the  order  of  103,  and 
turbulent  transport  is  unimportant  except  in  a small  initial 
region  after  the  nozzle  exit  where  the  turbulence  is  the  decay- 
ing turbulence  field  left  over  from  the  nozzle  wall  boundary 
layers  and  the  recirculation  regions  at  the  nozzle  exits. 
Chemical  laser  configurations  now  being  investigated  for  future 
systems  are  expected  to  operate  at  higher  cavity  pressures  and 
in  flow  regimes  where  turbulent  transport  processes  will  be 
more  important.  The  use  of  boundary  layer  trips  to  produce 
turbulent  boundary  layers  in  the  nozzles  has  led  to  some 
improvement  in  laser  performance,  but  some  of  the  observed 
effects  are  not  easily  explained.  In  some  cases  a decrease  in 
laser  performance  has  resulted  from  attempts  to  introduce  tur- 
bulence into  the  flow.  The  effects  of  various  turbulence 
parameters  on  the  mixing  and  reaction  in  the  chemical  laser 
are  quite  complex. 

The  flow  parameters  of  the  DF  chemical  laser  under  study 
have  been  shown  in  Figure  8.  The  cavity  pressure  in  the  actual 
system  is  10  torr,  and  the  flow  Reynolds  number  is  about 
5 x 10^.  The  effect  of  the  following  four  turbulence  param- 
eters on  the  formation  of  excited  DF  species  in  the  initial 
region  has  been  studied. 

1.  Initial  disturbance  amplitude  u ’ u ' j_n^t  with 
isotropic  distribution  u 'u ' = v»v»  = w»w»  . 

2.  Turbulence  distribution  with  constant  turbulence 
kinetic  energy  q^  = u ' u ' + v* v'  + wTwT  . 

3.  Turbulence  macroscale  A . 
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A.  Inclusion  or  neglect  of  the  mixedness  correlation 
a'B'/al  in  the  chemical  source  terms. 

Similar  calculations  have  also  been  performed  for  a cavity 
pressure  of  100  torr  with  the  same  flow  geometry  to  study  the 
effect  of  raising  the  operating  pressure  and  the  flow  Reynolds 
number.  The  results  for  IQ  torr  cavity  pressure,  the  normal 
operating  conditions,  are  discussed  first  followed  by  the 
results  for  100  torr  pressure. 

Laser  Cavity  Pressure  = 10  Torr 


The  behavior  of  the  maximum  value  of  the  u’u'  correlation 
at  different  downstream  axial  stations  is  shown  in  Figure  9- 
The  figure  shows  the  results  for  calculations  started  with  two 
different  initial  levels  of  turbulence,  u ' u ' = .00^  , . 0^  . 
The  turbulence  distribution  is  isotropic  and  the  macroscale  has 
the  normal  value  A = .05  (y  ^ - y 25)  » where  ^.75  and  y.25 
are  the  values  of  the  normal  coordinate  where  u - U2/U1  - U2  is 
.75  and  .25,  respectively.  The  turbulence  field  here  is  just 
the  rapidly  decaying  turbulence  left  over  from  the  nozzle  wall 
boundary  layers  and  the  dead  water  recirculating  regions  at  the 
nozzle  exit.  The  initial  turbulence  level  can  be  increased  by 
the  use  of  boundary  layer  trips  of  various  designs.  Due  to  the 
low  flow  Reynolds  number  at  p = 10  torr  , the  turbulence 
decays  very  rapidly,  and  a laminar  transport  analysis  would 
appear  to  be  justified  for  most  of  the  flowfield.  However,  as 
the  figure  shows,  in  this  initial  region  the  turbulence  level 
at  a point  is  strongly  dependent  on  the  amplitude  of  the  initial 
disturbance,  and  it  is  shown  later  that  there  are  significant 
differences  in,  for  example,  DF(2)  formation  as  a function  of 
the  initial  turbulence  amplitude. 

Figures  10  and  11  present  the  results  with  and  without  the 
inclusion  of  turbulent  scalar  correlations,  like  the  mixedness 
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Figure  10.  Comparison  of  DF(2)  production  with  and  without  the  inclusion  of  the  turbulent 
scalar  correlations  in  the  chemical  source  terms.  X = 1.25  cm 


Figure  11.  Comparison  of  DF(2)  production  with  and  without  the  inclusion  of  the  turbulent 
scalar  correlations  in  the  chemical  source  terms.  X = 2.5  cm 


* r 


. 


i 


I 


i 


J 

* 


lJ 


1 i 

Li 


correlation  a’B'/aB  , in  the  chemical  source  terms.  The 
concept  can  be  illustrated  by  the  following  simplified  problem. 
Consider  the  expression  for  the  rate  of  change  of  a species  a , 
due  to  chemical  reactions  in  a turbulent  flow.  Assume  a simple 
one-step  reaction  a + 3 -*■  y . Then, 


da  . 
dt 


-kaB 


da  . 
dt 


-kaB  - 


ka’B'  + ak • B ' + Bk’a'  + k’a'B’ 


Many  of  the  previous  studies  of  turbulent  reacting  flows 
neglect  the  group  of  terms  within  the  brackets  and  simply  use 


da 

dt 


-kaB 


for  the  chemical  source  term  calculation.  This  is  analogous  to 
the  source  term  for  laminar  flow,  and  we  designate  this  procedure 
as  the  "laminar  chemistry"  approximation.  Another  error  in  this 
approximation  is  the  evaluation  of  k ; usually  one  takes 
k = k(T)  which  is  only  true  for  laminar  flow.  In  the  "turbulent 
chemistry"  runs,  k is  evaluated  from  the  pdf  for  "typical  eddy" 
model  runs  and  from  second-order  accurate  expansions  of  the 
Arrhenius  rate  expression  for  "secondary"  model  runs,  and  all 
the  terms  in  the  brackets  are  included  in  the  evaluation  of  the 
chemical  source  term.  In  the  current  studies  on  the  DP  chemical 
laser  reported  here  the  "secondary"  model  has  been  used  for  the 
higher-order  scalar  correlations.  Further,  as  the  temperature 
changes  in  the  laser  flowfield  are  quite  small,  temperature 
fluctuation  (and  therefore,  kf  ) effects  have  been  neglected. 
Therefore,  in  the  current  studies,  the  "turbulent  chemistry" 
formulation  corresponds  to  the  use  of  the  expression 
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and  the  "laminar  chemistry"  approach  involves  the  neglect  of 
the  mixedness  correlation. 

The  computations  in  Figures  10  and  11  show  the  transverse 
profiles  of  DF(2)  at  distances  1.25  cm  and  2.5  cm  downstream  of 
the  nozzle  exit.  The  results  for  the  other  vibrational  levels 
are  similar.  The  initial  disturbance  was  isotropic  and  of 
amplitude  u • u * = .04  and  scale  A . The  results  show  that  for 
p = 10  torr  , the  inclusion  of  the  fluctuating  chemical  source 
terms  has  a negligible  effect  on  the  formation  of  DF(v).  The 
turbulence-chemistry  interaction  effects  are  very  small  for 
these  conditions.  It  is  shown  later  that  for  larger  scale 
disturbances,  there  are  some  significant  turbulence-chemistry 
interactions  even  at  p = 10  torr  . 

The  effect  of  various  initial  turbulence  levels  on  DF(2) 
formation  is  shown  in  Figures  12  and  13.  Figure  12  shows  the 
profiles  of  DF(2)  at  2.5  cm  axial  distance  downstream  of  the 
nozzle  exit  for  laminar  flow  and  u ’ u * = . 004  and  .04  . 

The  Integrated  DF(2)  formation  in  the  normal  direction  is 
shown  in  Figure  13  as  a function  of  the  axial  position.  As 
remarked  earlier,  due  to  the  rapid  decay  of  the  turbulence 
the  actual  turbulence  level  in  the  initial  region  is  strongly 
dependent  on  the  initial  disturbance  level  and  results  in 
significant  changes  in  the  formation  of  the  excited  species. 
Higher  initial  disturbances  predict  the  formation  of  larger 
amounts  of  DF(2)  and  thus  should  lead  to  improved  laser  perfor- 
mance. Deactivation  reactions,  of  course,  must  be  included  in 
the  analysis  to  correctly  predict  the  laser  performance. 

The  use  of  boundary  layer  trips  in  the  nozzles  would  intro- 
duce initial  disturbances  of  variable  amplitude,  isotropy,  and 
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Figure  12.  Effect  of  various  initial  turbulence  levels  on 
DF(2)  formation. 
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scale.  The  effects  of  isotropy  and  turbulent  macroscale  vari- 
ation have  been  studied,  and  are  discussed  next. 

The  effect  of  turbulence  isotropy  was  investigated  by 
the  following  three  runs.  The  turbulent  kinetic  energy 
was  kept  constant  at  0.12  and  the  scale  A was  used.  The 
three  distributions  used  are  tabulated  below,  and  correspond 


to  typical  limiting  distributions  that  are  observed  in  a variety 
of  atmospheric  turbulent  flows.  The  results  for  p = 10  torr 
are  shown  in  Figure  it  for  two  axial  stations.  There  are  small 
differences  very  close  to  the  nozzle  exit,  but  in  the  region  of 
interest  for  the  laser  cavity  (~  2 cm)  there  is  virtually  no 
effect  of  the  initial  isotropy  characteristics  of  the  distur- 
bances. 


The  current  version  of  the  reacting  shear  layer  (RSL)  pro- 
gram uses  a constant,  algebraically  specified  turbulence  macro- 
scale A across  the  entire  flow  region  in  the  normal  direction. 
The  scale  is  based  either  on  the  mean  axial  velocity  profile  or 

I 

on  the  profile  of  the  turbulence  kinetic  energy.  In  all  the 
laser  calculations  discussed  here,  the  scale  is  based  on  the 
mean  velocity  profile  and  A = 0.5  (y.75  - y.25)  * We  are 
currently  engaged  in  the  development  of  a transport  equation 
for  the  turbulence  macroscale  for  compressible  flows  using  the 
formulation  for  a two-point  correlation  function.  Calculations 
have  been  performed  to  show  the  effect  of  scale  variation 
ranging  from  0.7  A to  2.0  A . The  scale  modification  parameter. 
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i.e.,  0.7  or  2.0 , remains  constant  in  the  program  calculations 
although  in  the  actual  flow  if  a larger  scale  disturbance  were 
introduced  at  the  initial  station,  it  would  adjust  to  the 
local  flowfield  conditions  in  some  characteristic  relaxation 
time  of  the  order  of  A/q  . Figure  15  shows  the  effect  of 
doubling  the  scale  on  the  transverse  profile  of  DF(2)  at 
the  axial  station  X = 2 . 5 cm  . The  result  is  not  surprising 
as  the  larger  scale  will  result  in  increased  transverse  mixing 
and  a more  uniform  profile.  The  calculations  were  done  using 
the  turbulent  chemistry  formulation.  It  has  been  shown  before 
that  for  10  torr  runs  with  scale  A , the  turbulent  chemistry 
and  the  laminar  chemistry  results  are  virtually  identical. 

However,  when  the  scale  is  increased  to  2A,  there  are  signifi- 
cant differences.  Figure  16  shows  the  results  for  10  torr, 

2A  runs  using  laminar  and  turbulent  chemistry  at  X = 2.5  cm  . 
Figure  17  presents  the  integrated  DF(2)  mass  fraction  as  a 
function  of  the  axial  position  for  the  same  conditions.  The 
comparison  shows  substantial  differences  between  the  laminar 
and  turbulent  chemistry  formulations,  and  significant  errors 
can  result  from  the  neglect  of  the  turbulence-chemistry  inter- 
action. The  behavior  of  the  species  mixedness  correlation 
a’B'/ a 3 for  10  torr  and  scales  A and  2A  is  shown  in  Figure  18. 
In  the  turbulent  chemistry  runs  we  normally  start  with  an  initial 
value  of  the  correlation  of  -1  across  the  mixing  layer.  For 
scale  A , the  mixedness  rapidly  drops  to  negligible  values  in 
line  with  the  behavior  of  the  velocity  turbulence  correlation, 
and  this  is  the  reason  why  Figures  10  and  11  showed  virtually  no 
difference  between  the  laminar  and  turbulent  chemistry  calcu- 
lations. When  the  scale  is  increased  to  2A,  the  mixedness 
correlation  has  fairly  substantial  values  for  the  entire  flow- 
field  and  this  leads  to  significant  differences  between  the 
laminar  and  turbulent  chemistry  calculations. 

The  overall  effect  of  scale  changes  on  laser  performance 
can  be  estimated  by  comparing  the  integrated  DF(v)  mass  fractions 
across  the  normal  profile.  Figure  19  presents  the  results  for 
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Figure  17.  Comparison  of  laminar  and  turbulent  chemistry 
formulations  on  integrated  DF(2)  formation. 


Integ 


scale  parameters. 


the  four  different  scale  parameters.  The  results  show  that  in 
the  initial  region  of  the  flow  (up  to  X = 2 cm  ) there  are 
substantial  differences  in  total  DF(2)  formation,  with  a smaller 
scale  leading  to  increased  laser  performance.  Further  down- 
stream, the  scale  effect  is  smaller  in  magnitude,  and  a cross- 
over of  the  predicted  curves  is  seen  indicating  improved 
performance  for  the  larger  scales. 

These  results  for  the  10  torr  cavity  pressure  DF  laser 
show  the  significant  effects  of  the  two  turbulence  parameters — 
initial  disturbance  amplitude  and  the  turbulence  macroscale  — on 
the  performance  of  the  laser.  Different  combinations  of  tur- 
bulence amplitude  and  scale  could  lead  to  improvement  or 
degradation  of  the  laser  performance  and  may,  in  part,  explain 
the  discrepancies  in  experimental  studies.  Of  course,  as 
remarked  earlier,  our  studies  must  be  repeated  with  a more 
complete  set  of  reactions  including  laser  radiation  effects 
to  enable  direct  comparison  to  experiments  and  to  make  the 
results  applicable  to  the  actual  laser  system. 

The  Important  results  for  the  10  torr  cavity  pressure  DF 
laser  are  summarized  in  Figure  20.  The  integrated  DF(2)  mass 
fraction  at  X = 1.25  cm  is  plotted  as  a function  of  the  initial 
turbulent  kinetic  energy  q2  , for  a number  of  different  scale 
runs.  There  are  only  five  data  points  on  the  graph,  but  the 
behavior  of  DF(2)  formation  with  q^  can  be  reasonably  expected 
to  have  the  trends  indicated  in  the  figure.  If  the  nozzle  wall 
boundary  layers  are  tripped  by  gas  injection  through  the  nozzle 
walls,  increases  in  turbulence  may  be  accompanied  by  increases 
in  the  scale,  and  it  may  be  hypothesized  that  for  increased 
blowing  one  follows  the  path  indicated  on  the  figure.  There 
is  an  initial  improvement  in  DF(2)  formation  (and  laser  perfor- 
mance), but  further  increase  in  blowing  results  in  an  adverse 
effect  on  the  total  DF(2)  production  due  to  the  increased  tur- 
bulence scale.  The  above  is  merely  a hypothetical  example  of 
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Figure  20.  Schematic  diagram  showing  decrease  in  DF(2)  production  due  to  boundary  layer 
tripping  by  blowing. 
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the  kind  of  effects  that  may  be  encountered  in  attempts  to 
modify  the  turbulence  parameters  in  the  chemical  laser  flowfield. 
There  is  very  little  information  available  on  the  scales  and 
intensities  of  turbulence  produced  by  various  trip  mechanisms 
and  some  investigation  of  this  should  be  carried  out.  Our 
program  calculations  clearly  demonstrate  the  effects  of  these 
turbulence  parameters  on  the  production  of  vibrationally  excited 
DF  species  in  the  initial  region  of  the  chemical  laser. 

Laser  Cavity  Pressure  = 100  Torr 

The  pressure  recovery  of  supersonic  mixing  chemical  lasers 
could  be  improved  by  operating  at  higher  cavity  pressures. 

However,  collisional  deactivation  also  increases  at  higher 
pressures  and  more  rapid  mixing  and  reaction  is  necessary  to 
maintain  or  improve  the  laser  performance.  Turbulence  provides 
a convenient  mechanism  to  increase  the  rate  of  mixing.  At 
higher  cavity  pressure  the  flow  Reynolds  number  is  higher  and, 
therefore,  it  is  possible  to  operate  in  a regime  where  turbulent 
transport  processes  would  be  dominant.  The  reacting  shear 
layer  program  has  been  used  to  study  the  effects  of  raising  the 
cavity  pressure  to  100  torr. 

The  behavior  of  the  velocity  turbulence  correlations  is 
shown  in  Figure  21.  The  flow  Reynolds  number  is  of  the  order 
of  5 x 10^  , and  the  turbulence  approaches  an  equilibrium  value 
in  a distance  of  a few  centimeters.  The  initial  amplitude  of 
u ' u * does  not  affect  this  equilibrium  value  and,  as  will  be 
shown  later,  does  not  significantly  affect  the  total  DF(2) 
formation.  The  behavior  of  uTuTmax  for  p = 10  torr  is  also 
shown  on  the  figure  for  comparison. 

Figure  22  shows  the  comparison  of  "laminar  chemistry"  versus 
"turbulent  chemistry"  at  p = 100  torr  for  two  axial  stations. 

All  the  100  torr  calculations  use  the  normal  turbulence  macroscale 
A . The  computations  show  very  substantial  differences  in  the 
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Figure  22.  Comparison  of  DF(2)  production  with  and  without  the 
inclusion  of  the  turbulent  scalar  correlations  in 
the  chemical  source  terms. 


DF(2)  formation.  The  effect  of  the  species  fluctuations  is 
to  reduce  the  rate  of  the  reaction  below  that  based  on  the 
means  alone,  for  the  species  are  not  mixed  on  the  molecular 
scale  to  the  extent  predicted  by  the  time-average  concentrations. 
For  very  rapid  reactions  the  species  mixedness  correlation 
cH  fT» / oi  B approaches  -1,  so  that  the  rate  of  reaction  could  be 
orders  of  magnitude  smaller  than  that  predicted  by  kaB  . The 
effect  of  the  temperature  (or  reaction  rate)  fluctuation  terms 
is  more  complicated  for  it  depends  on  whether  they  are  positively 
or  negatively  correlated  with  the  various  species.  However,  the 
figure  clearly  demonstrates  that  significant  errors  can  result 
from  the  neglect  of  the  turbulence-chemistry  interaction.  The 
behavior  of  the  species  mixedness  correlation  for  the  100  torr 
runs  is  shown  in  Figure  23.  The  computations  normally  start 
with  an  initial  value  of  a * 6 ' / a H = -1  across  the  mixing 
layer.  Some  runs  were  started  with  the  correlation  having  the 
initial  value  of  zero,  and  as  shown  in  the  figure  the  same 
general  behavior  is  observed.  The  10  torr  cavity  pressure 
results  for  scale  A are  also  plotted  for  comparison.  For 
the  10  torr  case  the  correlation  rapidly  drops  to  negligible 
values.  For  p = 100  torr  , ot'B'/a  8 has  fairly  substantial 
values  for  the  entire  flowfield  and  approaches  -1  at  large 
downstream  distances.  This  behavior  is  consistent  with  the 
results  presented  earlier  (Ref.  1)  showing  the  variation  of 
a'  B ' / a 6 with  Damkohler  number.  The  mixedness  correlation  in 
a free  shear  layer  approaches  -1  at  large  Damkohler  numbers. 

The  effect  of  the  initial  turbulence  amplitude  on  DF(2) 
formation  is  shown  in  Figures  2k  and  25.  Figure  2k  shows  the 
profiles  of  DF(2)  at  2.5  cm  distance  downstream  of  the  nozzle 
exit  for  laminar  flow  and  u fu » = *00k  and  .04  . The 
profiles  are  somewhat  different  but  the  peak  DF(2)  mole  fractions 
are  nearly  the  same.  The  integrated  DF(2)  formation  as  a 
function  of  the  axial  position  is  presented  in  Figure  25,  and 
indicates  somewhat  larger  production  of  excited  species  for 


59 


Figure  23.  Behavior  of  the  mixedness  correlation  at  two  different  cavity  pressures 
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Figure  24.  Effect  of  various  initial  turbulence  levels  on  DF(2)  formation 
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Figure  25.  Effect  of  various  initial  turbulence  levels  on 
integrated  DF (2)  production. 


larger  Initial  disturbance  amplitudes.  The  effect  of  scale 
variations  was  not  investigated  for  the  100  torr  cavity  pressure. 

The  effect  of  initial  turbulence  isotropy  was  investigated 
for  the  100  torr  pressure  laser  in  the  same  manner  as  for  the 
10  torr  case,  and  similar  results  were  obtained.  There  are 
small  differences  in  the  DF(2)  predictions  very  close  to  the 
nozzle  exit,  but  further  downstream  there  are  only  negligible 
differences  among  the  cases  that  were  investigated. 

Figure  26  compares  the  predictions  for  DF(2)  formation  at 
X = 2.5  cm  for  p = 10  torr  and  p = 100  torr  . The  compu- 
tations have  the  same  initial  turbulence  amplitude  and  use  the 
full  "turbulent  chemistry"  formulation.  The  results  show  that 
significantly  larger  amounts  of  DF(2)  (and  other  excited  vibra- 
tion levels)  are  produced.  However,  as  the  complete  set  of 
deactivation  processes  is  not  included,  we  cannot  predict 
whether  or  not  this  translates  into  improved  laser  performance. 

Figure  27  shows  the  typical  axial  variation  of  DF(v)  and 
the  temperature  in  the  central  region  of  the  flow  for  100  torr 
cavity  pressure.  These  results  were  calculated  using  the 
turbulent  chemistry  formulation,  and  the  reaction  scheme  in 
Table  1.  The  results  were  somewhat  surprising  in  the  fact  that 
the  predictions  showed  DF(3)  and  DF(2)  larger  than  DF(1).  The 
runs  were  repeated  with  the  more  complete  reaction  scheme  of 
Table  2 and  virtually  the  same  results  were  obtained.  To 
further  verify  the  basic  correctness  of  the  RSL  program  and 
the  multi-step  reaction  procedure  being  used  in  the  program, 
the  calculations  were  also  repeated  using  an  eddy  viscosity, 
laminar  chemistry  code  available  at  A.R.A.P.  (Low  Altitude  Plume 
Program  or  LAPP),  and  the  same  trends  in  the  DF(v)  axial  profiles 
were  predicted.  Therefore,  the  Reacting  Shear  Layer  program 
calculations  with  the  use  of  the  reaction  schemes  outlined  in 
Tables  1 and  2 appear  to  be  correct.  Calculations  should  be 
performed  with  the  use  of  a more  complete  reaction  scheme  in 
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Figure  26.  Effect  of  cavity  pressure  on  DF(2)  production 


our  program  for  detailed  comparison  with  experimental  measure- 
ments in  chemical  laser  systems. 


These  calculations  for  DF  chemical  lasers  using  our 
complete  second-order  closure  program  have  clearly  demonstrated 
the  significant  effects  of  the  turbulence  characteristics  on 
the  performance  of  chemical  lasers  operating  in  the  transitional 
flow  regime.  Significant  turbulence-chemistry  interactions 
could  be  present  at  somewhat  higher  flow  Reynolds  numbers.  The 
above  calculations  have  to  be  repeated  with  the  "typical  eddy" 
model  to  check  the  sensitivity  of  the  results  to  the  pdf  models 
for  the  higher-order  scalar  correlations.  Further,  the  more 
complete  reaction  set  of  laser  kinetics  should  be  used  in  the 
program  to  compare  the  program  predictions  to  experimental 
measurements  and  to  predict  the  effect  of  turbulence  parameters 
on  actual  laser  systems. 


VI.  SUMMARY  AND  CONCLUSIONS 


A complete  second-order  closure  program  for  the  study 
of  turbulent  reacting  shear  flows  has  been  developed.  The 
program  Is  capable  of  handling  multi-species , multi-step 
chemical  reactions  and  includes  turbulence-chemistry  interaction 
effects  in  the  analysis. 

The  program  has  been  used  to  study  the  mixing  and  chemical 
reactions  in  the  initial  region  of  a supersonic  DF  chemical 
laser.  The  effects  of  a number  of  turbulence  variables  — 
intensity  q^  t scale  A , and  the  mixedness  correlation 
a'S'/a  B — on  the  formation  of  vibrationally  excited  DF(v) 
species  have  been  studied.  The  studies  at  low  pressure 
(10  torr)  demonstrate  the  need  for  a second-order  approach 
to  predict  the  details  of  the  transitional  turbulent  flow. 

There  are  very  significant  effects  of  the  initial  turbulence 
amplitude  and  turbulence  scale  on  the  formation  of  DF(v). 
Different  combinations  of  q2  and  A could  lead  to  improve- 
ment or  degradation  in  laser  performance.  At  higher  cavity 
pressures  (100  torr),  the  flowfield  has  a significant  level 
of  turbulence  and  there  are  Important  turbulence-chemistry 
interaction  effects  that  must  be  taken  into  account.  Signif- 
icant errors  can  result  from  the  neglect  of  the  mixedness 
correlations  in  the  chemical  source  terms  for  finite  rate 
chemistry  calculations. 

The  important  conclusions  of  the  study  for  the  MESA  IV 
flow  parameters  are  summarized  in  the  table  below: 
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Pressure 

Effect  on  DF(v)  Formation  of: 

q2init 

Turbulence 

Isotropy 

A 

oPT5-/  a 5 

10  torr 

Strong 

Weak 

Strong 

Weak 

100  torr 

Strong 

Weak 

- 

Strong 

These  studies  have  clearly  demonstrated  some  of  the  important 
effects  in  nonequilibrium  turbulent  flowfields  that  can  only  be 
analyzed  using  a second-order  closure  program.  The  A.R.A.P. 
Reacting  Shear  Layer  (RSL)  program  has  now  reached  the  stage  of 
development  wherein  it  can  be  used  for  the  analysis  of  complex 
turbulent  reacting  flowfields. 


VII.  NOMENCLATURE 
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A 

b 

sij 

h 
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\i±  or  u,  v,  w 
Vc 

x , X 

y 

y. 75>  y.25 
a,  6,  Y 

6 

A 

\ 

P 

Superscripts 

— <> 

9 

t 


turbulence  model  parameter 
turbulence  model  parameter 

slit  width,  also  turbulence  model  parameter 

matrix  tensor 

sensible  enthalpy 

thermal  conductivity 
reaction  rate 

half  width  of  velocity  profile 
half  width  of  temperature  profile 
pressure 

turbulence  kinetic  energy 
turbulence  model  parameter 
temperature 
velocity  components 
turbulence  model  parameter 
axial  coordinate 
normal  coordinate 

position  in  flow  where  u - u2/ui  - u2  * . 75 , .25 
species  mass  fractions 
T — T 

1 -‘■ambient 
turbulent  macroscale 

turbulent  microscale 

density 

denotes  time  average 

denotes  fluctuation  about  the  mean  value 
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TABLE  1,  DF  CHEMICAL  LASER  REACTION  SET 


REACTIONS 


Kp  = A*EXP(B/RT)/T**N 
(MOLES-ML-SEC  UNITS) 


M 

A 

N 

B 

1 

F 

+ 

D2 

= 

DF1 

+ 

D 

0.1929E  16 

0.9 

-2000 . 0 

2 

F 

+ 

D2 

= 

DF2 

+ 

D 

0.4582E  16 

0.9 

-2000.0 

1 

3 

F 

+ 

D2 

= 

DF3 

+ 

D 

0.7235E  16 

0.9 

-2000.0 

4 

F 

+ 

D2 

= 

DF4 

+ 

D 

0.5426E  16 

0.9 

-2000.0 

5 

DF1 

+ 

DFO 

= 

DFO 

+ 

DFO 

0.1025E  03 

-2.9 

3200.0 

6 

DF2 

+ 

DFO 

= 

DF1 

+ 

DFO 

0.1989E  03 

-2.9 

3200.0 

7 

DF3 

+ 

DFO 

= 

DF2 

+ 

DFO 

0.3015E  03 

-2.9 

3200.0 

8 

DF4 

+ 

DFO 

= 

DF3 

+ 

DFO 

0.4040E  03 

-2.9 

3200.0 

9 

DF1 

+ 

D2 

= 

DFO 

+ 

D2 

0.5607E  02 

-3.0 

0.0 

B 

i 10 

DF2 

+ 

D2 

= 

DF1 

+ 

D2 

0.1145E  03 

-3.0 

0.0 

*- 

( 

1 n 

DF3 

+ 

D2 

= 

DF2 

+ 

D2 

0.1688E  03 

-3.0 

0.0 

\ 

12 

DF4 

+ 

D2 

= 

DF3 

+ 

D2 

0.2231E  03 

-3.0 

0.0 

h 

13 

DF1 

+ 

F 

= 

DFO 

+ 

F 

0.2834E  16 

0.7 

-3600.0 

K . 

14 

DF2 

+ 

F 

= 

DF1 

+ 

F 

0.3316E  16 

0.7 

-3600.0 

i 

15 

DF3 

+ 

F 

= 

DF2 

+ 

F 

0.4462E  16 

0.7 

-3600.0 

16 

DF4 

+ 

F 

= 

DF3 

+ 

F 

0.8441E  16 

0.7 

-3600.0 
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TABLE  2.  DF  CHEMICAL  LASER  REACTION  SET 


REACTIONS 


1 

F 

+ 

D2 

— 

DF1 

+ 

D 

2 

F 

+ 

D2 

= 

DF2 

+ 

D 

3 

F 

+ 

D2 

= 

DF3 

+ 

D 

A 

F 

+ 

D2 

= 

DF4 

+ 

D 

5 

DF1 

+ 

M 

= 

DFO 

+ 

M 

6 

DF2 

+ 

M 

= 

DF1 

+ 

M 

7 

DF3 

+ 

M 

= 

DF2 

+ 

M 

8 

DF4 

+ 

M 

= 

DF3 

+ 

M 

WHERE 

M 

5 

> DF(v> 

0- 

-4 

9 

DF1 

+ 

D2 

= 

DFO 

+ 

D2 

10 

DF2 

+ 

D2 

= 

DF1 

+ 

D2 

11 

DF3 

+ 

D2 

= 

DF2 

+ 

D2 

12 

DF4 

+ 

D2 

= 

DF3 

+ 

D2 

13 

DF1 

+ 

F 

= 

DFO 

+ 

F 

14 

DF2 

m 

+ 

F 

= 

DF1 

+ 

F 

15 

DF3 

+ 

F 

= 

DF2 

+ 

F 

16 

DF4 

+ 

F 

= 

DF3 

+ 

F 

17 

DF1 

+ 

HE 

= 

DFO 

+ 

HE 

18 

DF2 

+ 

HE 

= 

DF1 

+ 

HE 

19 

DF3 

+ 

HE 

= 

DF2 

+ 

HE 

20 

DF4 

+ 

HE 

= 

DF3 

+ 

HE 

Kp  = A*EXP(B/RT)/T**N 


(MOLES-ML-SEC  UNITS) 

A 

N 

B 

0.1929E  16 

0.9 

-2000.0 

0.4582E  16 

0.9 

-2000.0 

0.7235E  16 

0.9 

-2000.0 

0.5426E  16 

0.9 

-2009.0 

0.1025E  03 

-2.9 

3200.0 

0.2050E  03 

-2.9 

3200.0 

0.3075E  03 

-2.9 

3200.0 

0.4100E  03 

-2.9 

3200.0 

0.5607E  02 

-3.0 

0.0 

0.1145E  03 

-3.0 

0.0 

0.1688E  03 

-3.0 

0.0 

0.2231E  03 

-3.0 

0.0 

0.2834E  16 

0.7 

-3600.0 

0.3316E  16 

0.7 

-3600.0 

0.4462E  16 

0.7 

-3600.0 

0.8441E  16 

0.7 

-3600.0 

0.1561E-05 

-5.0 

0.0 

0.3123E-05 

-5.0 

0.0 

0.4685E-05 

-5.0 

0.0 

0.6271E-05 

-5.0 

0.0 
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APPENDIX  A. 


GENERAL  CONSERVATION  EQUATIONS  FOR  A TURBULENT 
MULTICOMPONENT  REACTING  SYSTEM 


Continuity  Equation 

pt  + 

(pu\d  = 0 

(1) 

Momentum  Equation 

PUl  + 

Jt 

puS,* = -p,j + tj\* 

(2) 

where 

d dk  / , „ \ 

Tj  = g U(Ui,k  + UkJ 

+ 6ju*uk.  (3) 

J > K 

Energy  Equation 

pht  + pu\d  " Pt 
where 


<J>  = x 


d 

,m 


(A) 

(5) 

(6) 

(7) 

(8) 


H = -g 


dm 


kT,l  + ‘U’Wa.l 


(9) 


.m 


Species  Equation 


pc„  + pu  c 


a , d 


w + g^m  fpD  c - \ 
a & v a a ,1 J, 


m 


Equation  of  State 


w + G^ 
a a 


p = pftT  l ^ 


(10) 

(11) 


A 


A-l 


In  the  above-noted  equations,  the  symbols  have  their  usual 


meaning;  w^  is  the  mass  rate  of  formation  of  component  a per 


unit  volume;  h^  is  the  heat  of  formation  of  species  a ; D&  is 


the  diffusion  coefficient;  and  is  the  species  molecular 


weight . 

Writing  the  dependent  variables  in  these  equations  as  the 
sum  of  a mean  and  a fluctuating  part  (u  = u + u'  ; p = p + p'  ; 
etc.),  one  can  derive  equations  for  the  mean  quantities  and  the 


second-order  correlations  of  the  fluctuations. 


-Z, 


Defining  >yf=f.  + uf#  , the  equations  can  be  written 

T/  y 36 


as 


V P + pu££  + (p'u'  £)  ^ = 0 


(12) 


Z - 


P£>lL  + p'u'  0^  ^ + (p  'uj  )t  + 


p(u'£uj)  + u£(p'uj)  + (p'u'£uj) 


> 0- 


= "P»J  + 


(13) 


Z 7 


p£h  + p'u'  h t + ( p ' h ' )t  + 


p(u'  h')  + u (p'h')  + (p'u'  h') 


,*■ 


= p + U£P  + (u,AP;£)  + * - l ha«a 

" * rv 


(14) 


+ (p'u'  £)ca^  + (p'o^)t  + 


a 


p(u'£c')  + u£(p'c')  + (p'u'£c') 


= w + G 
a a 


(15) 


A-2 


A-5 


I i , 

iitikJ  l 


pap  (p’c^)  + (p'p'c^)t  + 


p(u'£p'c^)  + u£(p'p'c^)  + (p'p'u'£c^) 


J ,1 


- (p'u,£)  0(p'c’)  + (p'pMfic,  + (p'c')-^p 


,SL^  a 


a 


'a , SL 


'a /KK,£ 


a , a 


,z 


+(p’w^)  + (p’G^)  - 


pp(u;k}  + 2pu%(p,cc;) 


+ 2p(p'c^u|£)  + u££(p'p'c^)  + (p»p'c^u'£) 


(25) 


Detailed  expressions  for  terms  like  \p  , ip’  , G . G'  , etc, 

’ * ot  ’ a ’ 


can  be  obtained  from  their  definitions.  The  expressions  for  the 


reaction  source  terms  w and  w'  will  be  obtained  from  the 

a a 


modeling  of  the  chemical  reaction.  As  an  example,  if  one  uses 
a one-step  chemical  reaction  model 


} v 'M  — ► > v ' 'M 
ql  a s a a 

a=l  a=l 


(2  6) 


the  reaction  rate  term  can  be  written  as 


N / „ \v  ’ 
pc  \ a 


= W_(v”  - v:)  ATnexp  (-  J*  ~vr~ 


a a a 


a 


V «T/0  = 1 \ W 


(27) 


Expressions  for  wa  and  w^  can  be  obtained. 


A-7 


+ (u,£p')pc  0 + (u'£c')pp  + ( p ' p ' u 1 £ ) c . + ( p ' c ' u ' £ ) p „ 


